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Abstract 

We review some results obtained in the context of the Collaborative Research Center/Transregio 9. In particular 
we discuss three-loop corrections to the Higgs boson mass in the Minimal Supersymmetric Standard Model, higher 
order corrections to Higgs boson production, and the calculations of renormalization group functions and decoupling 
constants. 
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1. Introduction 

The discovery of a Higgs boson in July 2012 at the 
LHC has given a big boost to particle physics, both 
to theoretical developments and to experimental ana¬ 
lyses. On the experimental side an intensive study of the 
properties of the newly discovered particle has started. 
Among them are measurements of the mass, the spin, 
the couplings to fermions and bosons, the self-coupling 
and the decay width. At the same time the search for 
“New Physics”, i.e., phenomena that are not described 
by the Standard Model (SM) of particle physics, have 
been intensified. The experimental efforts are suppor¬ 
ted by theoretical studies which concentrate on the one 
hand on the development of new theories which can be 
tested experimentally. On the other hand, higher order 
quantum corrections are computed which are necessary 
to match the precision reached by the experimental col¬ 
laborations. In this review we discuss several calcu¬ 
lations in the context of the Higgs boson which have 
been performed within the Collaborative Research Cen¬ 
ter/Transregio 9 (CRC/TR 9). Some of the calculations 
are performed within the framework of the SM, others 
in extensions like the Minimal Supersymmetric Stand¬ 
ard Model (MSSM). 

Within the MSSM there are five physical Higgs bo¬ 
sons, two CP-even, one CP-odd and a charged one. A 


central feature of the MSSM is the prediction of the 
mass of the lightest CP-even Higgs boson. At lowest or¬ 
der it is bounded by the Z boson mass but higher order 
corrections, which are significant, can raise the value to 
above 130 GeV. Thus, there are parameter sets of the 
MSSM which are consistent with the Higgs boson mass 
value of about 125 GeV as measured at the LHC. In 
Section [2] we describe the code H3m which is the only 
publicly available program containing complete strong 
three-loop corrections. 

A crucial input in the experimental analyses con¬ 
cerned with Higgs boson properties are precise predic¬ 
tions of its production cross section. Section[3]discusses 
several calculations in this context. Among them are the 
next-to-next-to-leading order (NNLO) QCD corrections 
performed in full theory. Taking the limit of large top 
quark mass allows a quantitative check of the effective- 
theory result which is implemented in most of the com¬ 
puter codes. Furthermore, we summarize the first steps 
towards the third-order corrections which are needed for 
the phenomenological analyses since sizeable correc¬ 
tions are expected. Moreover, this calculation is quite 
challenging from the technical point of view. Section [3] 
also contains a summary of supersymmetric (SUSY) 
corrections to Higgs boson production in the gluon fu¬ 
sion channel, in particular, three-loop corrections to the 
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matching coefficient of the effective Higgs-gluon coup¬ 
ling. 

Finally, Section[4]is devoted to renormalization group 
functions and decoupling relations. A particular em¬ 
phasis is put on the recent calculation of the three- 
loop corrections to the beta functions of the SM coup¬ 
lings. Furthermore, we describe the decoupling proced¬ 
ure which is relevant when crossing flavour thresholds 
within the SM but also for the transition from the SM 
to the MSSM and from the MSSM to a Grand Unified 
Theory (GUT). As applications we discuss low-energy 
theorems which relate decoupling constants with effect¬ 
ive Higgs boson couplings, gauge coupling unification 
in the minimal SUSY SU(5) model, and the stability of 
the Higgs potential in the SM. 


2. Higgs boson mass in the MSSM 


M h , is bounded from above. At tree-level, the mass mat¬ 
rix of the neutral CP-even Higgs bosons h , H has the 
following form: 


M 


2 

H, tree 


sin 2/3 
--X 


/ M 2 cot/3 + M 2 tan/3 

l -m\-K 


-Ml- 

M 2 tan y 3 + M 2 A cot/3 


( 1 ) 


The diagonalization of A4^ tree gives the tree-level res¬ 
ults for Mh and M H , and leads to the well-known bound 
Mh < Mz which is approached in the limit tan/3 —> oo. 
Radiative corrections to the Higgs pole mass raise this 
bound substantially. The dominant radiative corrections 
are generated by the top quark and top squark loops that 
scale like ~ a t m 2 ~ ( m t is the top quark mass and 

yfa~ t is proportional to the top-Yukawa coupling) 16}® • 
Including higher order corrections, one obtains for 
the Higgs boson mass matrix 


The Higgs boson mass measurement by ATLAS 
M h = 125.36 ± 031(stat.) ± 0.1S(syst.) GeV ffl 
and CMS M h = 125.02^(^.)^;{4( l yy^.) GeV El 
already reached the accuracy level of a precision observ¬ 
able. All the other properties (couplings, spin and par¬ 
ity) of the new particle have been determined with sig¬ 
nificantly lower accuracy. They are consistent with the 
SM predictions based on a minimal Higgs sector. How¬ 
ever, also, various other beyond-the-SM (BSM) theories 
with a richer Higgs sector can accommodate the present 
experimental data | liEl- While within the SM the Higgs 
boson mass is a free parameter, in BSM theories it can 
often be predicted, providing an important test of the 
model. For example, the mass of the lightest Higgs bo¬ 
son within supersymmetric models is, beyond the tree- 
level approximation, a function of the top squark masses 
and mixing parameters. It grows logarithmically with 
the top squark masses and can be used to determine 
the supersymmetric mass scale, once the mixing para¬ 
meters are fixed. This approach has received consider¬ 
able attention recently EHU, partially because the dir¬ 
ect searches for supersymmetric particles at the LHC re¬ 
mained unsuccessful, indicating a possible lower bound 
for the supersymmetric mass scale in the TeV range. 

In the following, we concentrate on the precise pre¬ 
diction of the lightest Higgs boson mass within the 
MSSM. Compared to the SM, the MSSM Higgs sec¬ 
tor is described by two additional parameters, usually 
chosen to be the pseudo-scalar mass Ma and the ra¬ 
tio of the vacuum expectation values of the two Higgs 
doublets, tan/3 = V 2 M. The masses of the other Higgs 
bosons are then fixed by supersymmetric constraints. In 
particular, the mass of the light CP-even Higgs boson, 


K = M^ trcc - ( V ¥*). (2) 

which again gives the physical Higgs boson masses 
upon diagonalization. The renormalized quantities , 
and 11^02 are obtained from the self energies of 
the fields 0i, 02 and A, evaluated at zero external mo¬ 
mentum, as well as from tadpole contributions of 0 i and 
02- Here 0i, 02 and A denote the CP-even and CP-odd 
neutral components of the two Higgs doublets. 

The one-loop corrections to the Higgs pole mass are 
known without any approximations EHID. The bulk of 
the numerical effects can be obtained in the so-called 
effective-potential approach, for which the external mo¬ 
mentum of the Higgs propagator is set to zero. Most of 
the relevant two-loop corrections have been evaluated in 
this approach (for reviews, see e.g. Refs. Ifl2lfl3ln . Ex¬ 
act calculations at two-loop order JT41 (15J showed ex- 
plicitely, that momentum-dependent contributions can 
reach the current experimental accuracy only for very 
heavy supersymmetric particles with masses above few 
TeV. In addition, two-loop corrections including even 
CP-violating couplings and improvements from renor¬ 
malization group considerations have been computed in 
Refs. If72l fT3l fT6l. In particular CP-violating phases 
can lead to a shift of a few GeV in Mh, see, e.g., 
Refs. (I3II8). In Ref. fl9l a large class of sub-dominant 
two-loop corrections to the lightest Higgs boson mass 
have been considered. 

The first complete three-loop calculation of the lead¬ 
ing quartic top quark mass terms within SUSY QCD has 
been performed in Refs. l20l [2H (see also Ref. l22l for 
a recent review). Given the many different mass para¬ 
meters entering the formula for the Higgs boson mass in 
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Eq. ^ an exact calculation at the three-loop level is cur¬ 
rently not feasible. However, it is possible to apply ex¬ 
pansion techniques 1231 for various limits which allow 
to cover a large part of the supersymmetric parameter 
space. For the case of the Higgs mass corrections, the 
occurring Feynman integrals can be reduced to three- 
loop tadpole topologies that can be handled with the 
program MATAD EQ Concerning renormalization, it 
is well known that the perturbative series can exhibit 
a bad convergence behaviour in case it is parametrized 
in terms of the on-shell top quark mass. This feature 
is due to intrinsically large contributions related to the 
infra-red behaviour of the theory. Thus, in Ref. lf2~T1l the 
results for the Higgs boson mass are expressed in terms 
of the top quark mass renormalized in the DR scheme. 
In addition, in order to avoid unnatural large radiative 
corrections for scenarios with heavy gluinos, a modified 
non-minimal renormalization scheme for the top squark 
masses has also been introduced in Ref. ED. The addi¬ 
tional finite shifts of top squark masses are chosen such 
that they cancel the power-like behaviour of the gluino 
contributions. 

For heavy supersymmetric particles (with a typical 
mass scale M S usyX the radiative corrections to the 
Higgs boson mass contain large logarithms of the form 
ln(m f /MsusY)- They have to be resummed in order 
to extend the validity range of the perturbative expan¬ 
sion up to large Msusy values. The dominant contri¬ 
butions to the leading (FF) and next-to-leading logar¬ 
ithmic (NFF) terms up to the fourth loop order have 
been obtained in Ref. f25l . Very recently, the gener¬ 
alization of the FF and NFF approximation to the sev¬ 
enth loop order has been derived n Furthermore, in 
Ref. 0 the recent calculations of the three-loop beta- 
functions for the SM coupling^] and the two-loop cor¬ 
rections to the Higgs boson mass in the SM [26l have 
been used to derive the (presumably) dominant next-to- 
next-to-leading logarithmic (NNFF) corrections at the 
four-loop order. Finally, the resummation of the large 
logarithms contained in the running of the top quark 
mass at three loops within SUSY QCD has been per¬ 
formed in Ref. l27ll and implemented in the code H3m. 

Usually, the resummation is achieved with the help 
of Renormalization Group Equations (RGEs). For ex¬ 
ample, the resummation performed in Ref. [27]] amounts 
to the following steps: the running top quark mass is 
determined in the SM with the highest available preci¬ 
sion EEH33 from the pole mass. Then, the running 
mass is evolved up to the scale where supersymmetric 


! For details see Section^ 


particles become active (SUSY scale) using the RGEs of 
the SM. Afterwards, the running top quark mass in the 
SM is converted to its value in the MSSM. In this step, 
threshold corrections at the SUSY scale are required. In 
the last step, the running top quark mass in the MSSM 
is evolved to the desired energy scale with the help of 
MSSM RGEs. As the RGEs in the SM and the MSSM 
are known to three-loop order, the threshold corrections 
are required at the two-loop order. The explicit calcu¬ 
lations of the three-loop RGEs and the two-loop super- 
symmetric threshold corrections are presented in detail 
in Section [4] 

There are by now several computer programs pub¬ 
licly available which include most of of the higher or¬ 
der corrections to the lightest Higgs boson mass in the 
MSSM. FeynHiggs has been available already since 
1998 M ED and has been continuously improved 
since then. In particular, its last version contains all 
numerically important two-loop corrections, as well 
as the resummation of the FF and NFF computed in 
Ref. (4). A second program, CPSuperH (32], is based 
on renormalization group improved calculations and al¬ 
lows for explicit CP violation. The third program, 
H3mED, contains all currently available three-loop res¬ 
ults together with the resummation of the large log¬ 
arithms derived in ED- In the latest version of H3m 
also the Mathematica package SLAM [33] has been 
implemented. SLAM provides an interface for calling 
and reading output from SUSY spectrum generators 
fully automatically, using input parameters specified in 
the Supersymmetry Fes Houches Accord (SFHA) f34l . 
Furthermore, the Higgs boson masses are also calcu¬ 
lated by the SUSY spectrum generators Sof tSusy EU, 
SPheno [36], and SuSpect [37] using DR paramet¬ 
ers and two-loop RGEs. For low supersymmetric mass 
scales below 1 TeV, the predictions for Mh of all codes 
are in quite good agreement with differences within 
about 1 GeV. This level of agreement is in general suf¬ 
ficient due to a sizeable parametric uncertainty for the 
prediction of Mh , that is mainly induced by the experi¬ 
mental uncertainty in the measurement of the top quark 
mass. However, for large SUSY scales in the multi- 
TeV range, the differences between the predictions of 
the various codes becomes substantial. This behaviour 
can be explained by the increase of the radiative correc¬ 
tions with the SUSY scale, and by the fact that differ¬ 
ent orders in perturbation theory are implemented in the 
various programs. 

For illustration we show in Fig. [T] a comparison of 
the predictions for the Higgs boson mass computed 
with H3m at two and three loops (long-dashed and solid 
lines), and FeynHiggs at two loops with and without 
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Figure 1: Comparison of the predictions for the lightest Higgs boson mass provided by H3m and FeynHiggs. The long-dashed and solid lines depict 
the predictions of H3m at two- and three-loop orders. The dashed line shows the two-loop result obtained with FeynHiggs including resummation of 
LL and NLL. The dotted line corresponds to the FeynHiggs result without resummation. The error bands account for the theoretical uncertainties 
of SUSY-QCD results induced by the variation of the renormalization scale from m t ! 2 to Msusy (see text for more details). 


resummation of LL and NLL (dashed and dotted lines). 
The experimentally measured value for the Higgs boson 
mass is depicted by the horizontal black band. 

For the supersymmetric input parameters we choose 
a constrained MSSM scenario with the following val¬ 
ues: Aq = 0, tan/3 = 10, /i > 0 and mo = mi /2 = 
200 GeV,..., 5 TeV where the definition M S usy = 
^/M stl M st2 has been used. M stl and M st2 denote the stop 
quark masses. For the results obtained with FeynHiggs 
and for the central values of the H3m predictions we fix 
the renormalization scale to the top quark pole mass. 
The discrepancy between the two-loop predictions of 
the two codes can be traced back to differences in the 
renormalization schemes and resummation procedures. 
For the MSSM scenario presented here the FeynHiggs 
prediction lies quite close to the central curve of H3m 
results at three loops. Note, however, that there are also 


other choices of MSSM parameters where bigger differ¬ 
ences between the resummed two-loop predictions of 
FeynHiggs and H3m are observed. 

The uncertainty bands around the H3m predictions 
have been obtained by varying the renormalization 
scale between m t /2 and M S usy using for simplicity the 
SUSY-QCD corrections up to three-loop order in the on- 
shell scheme as derived in Ref. (20). We are aware that 
there is an inconsistency with the central values (long- 
dashed and solid line) which are obtained in the mod¬ 
ified DR scheme ED Nevertheless, the bands should 
reflect the realistic size and variation of uncertainties. 
Indeed it is nice to see that the three-loop band lies al¬ 
most completely inside the two-loop band. 

As can be seen from the figure the radiative cor¬ 
rections increase with the SUSY mass scale and can 
amount even at the three-loop level up to few GeV. It 
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is important to mention that for heavy SUSY scales, the 
genuine three-loop radiative corrections are few times 
bigger than the parametric uncertainty induced by the 
uncertainty in the top quark mass measurement (estim¬ 
ated to be of around 1 GeV [211) and about an order of 
magnitude bigger than the current experimental uncer¬ 
tainty on M h . Furthermore, guided by the size of the¬ 
oretical uncertainties at three loops, we conclude that to 
cope with the experimental precision on M/*, for heavy 
SUSY particles even the four-loop contributions are re¬ 
quired. 

Finally, the determination of the lightest Higgs boson 
mass including three-loop corrections relaxes the lower 
bound on the SUSY mass scale to about 4 TeV, greatly 
improving the prospects for supersymmetry discovery 
at the upcoming run of the LHC. Once again, this nu¬ 
merical analysis highlights the importance of improved 
theoretical calculations of M h to refine the implications 
of the Higgs boson discovery for constraining the super- 
symmetric models. 


3. Higgs boson production in the SM and MSSM 


A crucial input for the discovery of a Higgs boson 
at the LHC was the precise prediction of the cross sec¬ 
tion. In fact, there are several production mechanisms 
which are nicely summarized in the reviews (38lf40ll . In 
this section we concentrate on the gluon fusion channel 
which gives numerically the largest contribution. 

A promising approach to compute higher order cor¬ 
rections to Higgs boson production is based on the Lag- 
range density 

Xetf =- C\-G liV G flv , (3) 

v 4 


which describes the effective coupling of a Higgs bo¬ 
son to up to four gluons. In Eq. 0 G^ y is the 
gluon field strength tensor and C\ is the coupling (or 
matching coefficient) resulting from integrating out the 
heavy degrees of freedom from the underlying theory. 
Within QCD, C\ only depends on the top quark mass 
via In (p 2 /mf) where p is the renormalization scale. 
It has been computed up to three loops (NNLO) in 
Refs. (41-143 ] and the four-loop calculation has been 
performed in (44l [45). Using renormalization group 
techniques even the five-loop expression could be de¬ 
rived in Refs. |44|[45]], which, however, depends on the 
unknown coefficient of the fermionic part of the five- 
loop beta function. 

In the MSSM, C\ becomes a complicated function of 
all heavy mass scales and p, see Subsection 33 for more 
details. 


In Refs. [46-49] it has been demonstrated that the 
NNLO prediction of cr(pp —> H + X) within the 
effective-theory approach of Eq. 0 approximates the 
exact SM result with an accuracy below 1%, in particu¬ 
lar for Higgs boson masses around 125 GeV. This issue 
is discussed in Subsection 13. 11 Numerical NLO calcu¬ 
lations (50l suggest a similar behaviour in the MSSM. 

Although the Higgs production cross section is 
known to NNLO, the contribution from unknown higher 
orders is estimated to be of the order of 10% which 
asks for a N 3 LO calculation. We summarize the cur¬ 


we briefly discuss the production of a light Higgs boson 
within the MSSM. 

3.1. Top quark mass dependent results up to NNLO 

LO contributions to cr(pp —> H + X) have already 
been computed end of the seventies in Refs. EEHS) and 
also the NLO QCD corrections are available since al¬ 
most 20 years [[55] [56| taking into account the exact de¬ 
pendence on the top quark mass (see also Ref. 1571 for 
analytic results of the virtual corrections). NLO elec- 
troweak corrections have been computed in Ref. (58l 
and mixed QCD-electroweak corrections are considered 
in [59|. 

At LHC energies the NLO QCD corrections amount 
to 80-100% of the LO contributions which makes it 
mandatory to compute higher order perturbative correc¬ 
tions. Beginning of the century three groups have inde¬ 
pendently evaluated the NNLO corrections [60-63] in 
the limit of infinitely heavy top quark using the effective 
Lagrange density of Eq. 0. NNLO corrections to the 
production of a pseudo-scalar Higgs boson have been 
computed in Refs. (641165) . At NLO this approximation 
works surprisingly well, leading to deviations from the 
exact result that are less than 2% for M h < 2 M t (see, 
e.g., Ref. (66l ). However, one has to keep in mind that 
next to Mh and M t the partonic cross section also de¬ 
pends on the partonic center-of-mass energy s which, 
in principle, reaches up to the beam energy of the LHC 
questioning the validity of the assumption M t —> oo. 
Thus, it is necessary to investigate the validity of the 
large top quark mass approximation at NNLO. 

A further reason why the effective theory approach 
needs to be checked at NNLO is the fact that sev¬ 
eral improvements over the fixed-order calculation have 
been constructed. Among them is the soft-gluon re¬ 
summation to next-to-next-to-leading (671168) (see also 
Ref. (69) ) and next-to-next-to-next-to-leading [70-731 
logarithmic orders and the identification (and resumma¬ 
tion) of certain n 1 terms GH which significantly im¬ 
proves the perturbative series. 


rent status in Subsection |3.2| Finally, in Subsection 3.3 
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At NNLO an exact calculation, as it has been per¬ 
formed at NLO, is currently out of range and thus ap¬ 
proximation methods have to be used in order to estim¬ 
ate the effect of a finite top quark mass. In the analyses 
performed in Refs. ESm [75123 there are actually 
two ingredients which allow the systematic reconstruc¬ 
tion of an approximation for the partonic cross section: 

• It is suggestive to evaluate the Feynman diagrams 
contributing to Higgs production in the full the¬ 
ory applying an asymptotic expansion for large top 
quark mass. In the approach of Refs. [47, 49], 
where the imaginary parts of forward-scattering 
amplitudes have been considered, this requires the 
computation of four-loop Feynman diagrams as 
shown in Fig. [2] for the gluon- and quark-induced 
channels. In practice, one- or two-loop vacuum in¬ 
tegrals with mass scale M t and one- or two-loop 
box diagrams with massless external particles and 
a massive Higgs boson in forward-scattering kin¬ 
ematics have to be computed. As a result one ob¬ 
tains the partonic cross section in an expansion in 
the inverse top quark mass which should provide 
a good approximation for M 2 » M 2 , s. Four and 
five expansion terms in 1 /M 2 have been computed 
for a scalar and pseudo-scalar Higgs boson, re¬ 
spectively. 

• The second ingredient is the leading high-energy 
behaviour for the partonic cross section obtained in 
Refs. 114811751 and 1781 for the scalar and pseudo¬ 
scalar Higgs boson, respectively. At NLO it is 
given by a constant, at NNLO, however, the lead¬ 
ing term is proportional to ln(x) (with x = M^/s) 
and the constant is not known. 

In Refs. 175] [78] the high-energy behaviour has been 
matched to the infinite-top quark mass result. A system¬ 
atic study taking into account higher order 1 /M t results 
has been performed in Refs. [46-49]. To illustrate the 
procedure we show in Fig.[3|the partonic cross section at 
NLO as a function of x = M^/s for M h = 130 GeV [471. 
The quantity A^ is defined via l47l 

frgg-H+x = Alo (a£» + A® + ...) , (4) 

where Alo collects various constants and the exact top 
quark mass dependence. Lines with longer dashes in¬ 
clude higher order terms in the 1 /M t expansion which 
only converges up to the threshold at x « 0.14. Never¬ 
theless, the approximation constructed in Ref. 147 1 (dot¬ 
ted line) agrees well with the exact result (solid line, 
obtained from HIGLU f79l ) and leads to a negligible 



Figure 2: Sample forward scattering diagrams whose cuts correspond 
to the LO, NLO and NNLO corrections to gg —> H + X, qg —> H + X 
and qq —> H + X. Dashed, curly and thick (thin) solid lines represent 
Higgs bosons, gluons and top (light) quarks, respectively. 



Figure 3: Partonic NLO cross sections for the gg channel as function 
of x = for Mh = 130 GeV. The exact and approximated results 
are shown as solid and dotted lines and the dashed lines correspond 
to expansions in l/M t (long-dashed curves contain more expansion 
terms). 


deviation for the hadronic cross section. Indeed, even 
the difference in the hadronic cross section computed 
from the approximated result and the M t —> oo expres¬ 
sion (short-dashed straight line for x —> 0 in Fig. [3]) 
is below 1% at NLO for a scalar Higgs boson with 
Mh - 125 GeV and of the order of 6% for a pseudo¬ 
scalar Higgs boson with mass Ma = 300 GeV. 

A similar behaviour is observed at NNLO: the dif¬ 
ference between the effective-theory result and the had¬ 
ronic cross section computed for finite top quark mass 
is below 1% for Mh = 125 GeV and can amount to 
about 10% for a pseudo-scalar Higgs boson with mass 
M a = 300 GeV. Thus, for a scalar Higgs boson as ob¬ 
served at the LHC finite top quark mass effects are negli¬ 
gible beyond Born approximation for the inclusive cross 
section. 
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(a) (b) (c) 



(d) (e) (f) 


Figure 4: Sample Feynman diagrams in forward scattering kinematics 
contributing to the N 3 LO cross section for Higgs boson production 
in gluon fusion. All cuts through the Higgs boson line have to be 
considered. Dashed and curly lines correspond to Higgs bosons and 
gluons, respectively. The black blob indicates the effective Higgs- 
gluon vertex according to Eq. <[3j. 


3.2. Status at NNNLO and beyond 

The various contributions which have to be con¬ 
sidered for cr(pp —> H + X) at N 3 LO are shown in Fig. [ 4 ] 
where amplitudes with forward-scattering kinematics 
are shown. From these diagrams all cuts through the 
Higgs boson lines have to be computed. They include 
three-loop virtual corrections [see (a) and (b)], two-loop 
virtual corrections in association with a real emission 
of a parton [see (c)], squared contribution of the one- 
loop real-virtual corrections [see (d)], the one-loop vir¬ 
tual contribution with two real emissions [see (e)], and 
the contribution with three real emissions [see (f)]. In 
addition collinear counterterms for the parton distribu¬ 
tion functions have to be taken into account in order 
to arrive at an infrared finite quantity. Different groups 
have provided building blocks for the N 3 LO Higgs bo¬ 
son production cross section which will be briefly sum¬ 
marized in the following. 

• In Ref. m the four-loop corrections to the match¬ 
ing coefficient C\ of the effective Lagrangian 0 
have been constructed from the three-loop decoup¬ 
ling constant for the strong coupling constant with 
the help of renormalization group methods and 
a low-energy theorem (see also Section |4.8| ). In 
Refs. [44j 45.1 the result has been confirmed by 
an explicit calculation of the four-loop decoupling 
constant. 

• The three-loop corrections to the massless Higgs- 
gluon [cf. Fig. [JJa) and (b)] (and photon-quark) 
form factor have been obtained by two independent 
calculations 18011811 (see also Ref. (82); fermionic 
contributions to the photon-quark form factor have 
already been obtained in Ref. (83)). 


• The 0(e) contributions to the NNLO master integ¬ 
rals has been computed in Refs. [49,.84]|. 

• Results for the LO, NLO and NNLO partonic cross 
sections expanded up to order 6 3 , e 2 and e 1 , re¬ 
spectively, have been published in Refs. [85, 86). 

• All contributions from convolutions of partonic 
cross sections with splitting functions, which are 
needed for the complete N 3 LO calculation, are 
provided in Refs. (55,[57). The results of (85| have 
been confirmed in Ref. (86 1 . 

• The full scale dependence of the N 3 LO expression 
has been constructed in Ref. [86]. 

• The single-soft current to two-loop order has been 
computed in Refs. [88 , 89j] which is an important 
ingredient to the two-loop corrections with one ad¬ 
ditional real radiation [cf. Fig. 4jc)]. The latter 
have been computed in Refs. (901 (9T) . 

• The single-real radiation contribution which ori¬ 
ginates from the square of one-loop amplitudes 
[see Fig.|4jd)] has been computed to all orders in e 
in Refs. El EH. 

• The soft limit of the phase space integrals for 
Higgs boson production in association with two 
soft partons [cf. Fig. |4je)] were computed in 
Refs. 194, [95) , in the latter reference even to all 
orders in e. 

• The triple-real contribution to the gluon-induced 
partonic cross section [cf. Fig. 00] has been 
considered in Ref. (%1 in the soft limit, i.e. for 
y - 1 ~M 2 ls —> 0. In particular, a method has been 
developed which allows the expansion around the 
soft limit. Two expansion terms in y are provided. 

• Three-loop ultraviolet counterterms are needed for 
a s [97,;98] and the effective operator (99j. 

• Two leading terms in the threshold expansion for 
the complete N 3 LO total Higgs production cross 
section through gluon fusion has been presented in 
Refs. (94 i 100 HOD . In these references, for the 
first time, a complete third-order cross section has 
been constructed (although only for y —> 0) which 
constitutes an important step. For physical applic¬ 
ations probably more terms in the threshold expan¬ 
sion are necessary II100) . 

• A further activity concerns the development of sys¬ 
tematic approaches to compute the master integrals 
for cr(pp —» H + X ), see, e.g., Refs. [96, 1102) . 
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• Several groups have constructed approximate 
N 3 LO results for the total cross section taking into 
account information from soft-gluon approxima¬ 
tion and the high-energy limit l68l[7011721173111031 - 
105 1. Differences in the numerical results can 
partly be traced back to different procedures used 
for the resummation of higher order logarithmic 
contributions. 

Several results also apply to third order corrections to 
the Drell-Yan process; see, e.g., Ref. am 

3.3. Higgs boson production in the MSSM 

A Higgs boson with a mass of about 125 GeV is both 
consistent with the SM taking into account the avail¬ 
able precision data (see, e.g., Ref. il07l ) but also with 
supersymmetric theories, in particular the MSSM (see 
Section [2]). Thus, it is important to investigate also the 
quantum corrections in this theory. In the recent years 
several groups have provided significant contributions 
in this respect, mainly in the effective-theory framework 
which requires the computation of loop corrections to 
the matching coefficient C\ and hence only vacuum in¬ 
tegrals have to be considered. Important two-loop con¬ 
tributions have been calculated in Refs. I50lfl08im2l . 

At NLO there are also considerations in the full the¬ 
ory. In Ref. lf50l the total cross section has been com¬ 
puted numerically considering both the top and bottom 
sector and building blocks for a (semi) analytic full- 
theory calculation have been provided in Refs. 1 1131 , 
114 1; a complete calculation along these lines is still 
missing]^] 

Recently the effective-theory NLO corrections have 
been implemented in the publicly available computer 
code SusHi rim At NNLO the rough approxima¬ 
tion of Ref. [ 1117 1 has been implemented, i.e., the genu¬ 
ine supersymmetric corrections to C\ have been set to 
zero at three loops. A comprehensive summary of all 
available contributions towards a precise prediction of 
the Higgs boson production cross section is provided in 
Ref. [118’] which also contains a detailed discussion of 
the theoretical uncertainties. 

In the remaining part of this section we will discuss 
the NNLO corrections to C\ which have been computed 
within the CRC/TR 9. Sample diagrams contributing 
to Ci at one, two and three loops are shown in Fig. [5] 
The symbols t, t t , g, g , h and s denote top quarks, top 
squarks, gluons, gluinos, Higgs bosons and s scalars, 


2 See Ref. (TB) for preliminary results. 


respectively. The latter are auxiliary particles intro¬ 
duced to implement regularization by Dimensional Re¬ 
duction 11191 which respects supersymmetry. For the 
computation of Ci it is possible to expand the Feynman 
integrals in the external gluon momenta which leads to 
vacuum integrals. In contrast to the SM, in the MSSM 
many different mass scales are present which increases 
significantly the complexity of the calculation. In fact, 
the currently available tools do not allow for an exact 
calculation and one has to rely on approximation meth¬ 
ods. In Refs. ifHoirm relations between the various 
masses have been assumed such that phenomenologic¬ 
ally interesting scenarios can be studied. This includes 
both strong hierarchies among masses but also expan¬ 
sions in the mass differences. Note that in the latter 
case the expansion series can be written down in dif¬ 
ferent but equivalent ways. For example, for x ^ 1 
(where x stands for a ratio of masses) one can choose 
1 - x, 1 - 1/x, 1 - x 2 or 1 - 1/x 2 as expansion para¬ 
meter which all lead to the same result after considering 
the expansion up to a fixed order. However, in prac¬ 
tice different choices lead to different numerical results 
and also show different convergence properties (see also 
Ref. Il22l ). In the calculation of Refs. [ 1120 , 12111 up 
to five different masses occur which leads (for the hier¬ 
archies considered in Ref. ESP to up to 48 different 
possible representations. We have implemented soph¬ 
isticated expansion schemes with the purpose to select 
the representation with largest radius of convergence 
providing at the same time reliable error estimates for 
each point in the parameter space (see Refs. [ 120 , 112311 
for a detailed discussion of the algorithm). In this way 
three-loop corrections to C\ have been evaluated for the 
top and bottom sector, neglecting, however, the bottom 
Yukawa coupling. 

In the MSSM the coupling of the light CP-even Higgs 
boson to bottom quarks is proportional to tan/3 whereas 
the top quark-Higgs coupling is proportional to 1 / tan j3. 
Thus, for large values of tan/3 (in practice this means 
tan/3 >10; see, e.g., Ref. ill2l ) the bottom sector caij^] 
contribute significantly to Higgs boson production. The 
computation of the Feynman integrals with internal bot¬ 
tom quark loops are more challenging since it is not 
possible to apply an effective theory approach as for 
the top sector. Indeed, up to now only NLO correc¬ 
tions are available, both for the SM [56) and for SUSY 
PCD 15Q1ITTT1ITT21 . 

At this point we want to mention some field- 
theoretical issues connected to the s scalar which have 


3 Whether it indeed leads to large corrections depends on the values 
of the other parameters like the mixing angles in the Higgs sector. 
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Figure 5: Feynman diagrams contributing to C\ within the MSSM. 


been addressed in the course of the calculation per¬ 
formed in Ref. I1120M . In fact, besides a mass term for the 
s scalar also a coupling of two s scalars to a Higgs bo¬ 
son, which emerges through radiative corrections, has 
to be added to the Lagrange density. Its non-standard 
part then reads 

where s 0,a denotes the bare s scalar field and the dimen¬ 
sionful quantity A® mediates the coupling of the Higgs 
boson to s scalars. It is convenient to renormalize the £ 
scalar mass on-shell requiring M ®' s = 0, and A® via the 
following condition 

(A°) 2 = 6A 2 e + A 2 e , (6) 

where SAl is fixed via the condition that the renormal¬ 
ized coupling of the £ scalars to Higgs bosons is zero. 
Analytic results for the corresponding counterterms, 
which are needed up to the two-loop level in the case 
of A e , can be found in Ref. 1 1201 . 

In a first step a simplified scenario with degenerate 
supersymmetric masses has been considered in E1. 
In Ref. 111201 the results have been generalized by con¬ 


sidering various hierarchies of the involved supersym¬ 
metric particle masses. Furthermore, details on the 
renormalization procedure and the treatment of evanes¬ 
cent couplings have been discussed. The results of 1120] 
have been cross-checked in Ref. ED where a low- 
energy theorem has been used in order to obtain C\ from 
the decoupling constant of a s . 

In Fig. [6] we present results for the total cross section 
which is computed according to the formula 


cr(PP —> h + X) - (l + d EW ) X 


SQCDz x 

°v oo 

SQCDz x 

- cr t v 0 Us) 


NLO 


NLO 


+crf QCD 0u s , J u /! ) 


NNLO 


(7) 


where cr^ QCD refers to the NLO cross section includ¬ 
ing all top and bottom effects l lllll II1211 . After sub¬ 
tracting the top quark/top squark contributions with the 
help of ct^ qcd (ju^)| NL o we can add the result from the 
top quark/top squark up to NNLO. Note that cr t also 
contains numerically small contributions from a non¬ 
vanishing Higgs-bottom squark coupling whereas the 
bottom Yukawa coupling is set to zero, see Ref. jl2Qj 
for details. Finally, electroweak effects [58] are taken 
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into account in a multiplicative way. Note that they are 
only available in the SM and potential large MSSM ef¬ 
fects are neglected in Eq. 0- 

In Fig. [6] we discuss numerical effects of the in¬ 
dividual terms in Eq. (|7j) using the m™ ax scenario of 
Ref. (TSl as a basis. We apply slight modifications 
which lead to the following parameters (see Ref. IH201 
for explanations of the parameters) 


A b =A T = 2469.48 GeV, 
Mi = 5 s 2 w K3c 2 w )M 2 , 
M 3 = 800 GeV, 
yu susy = 200 GeV, 
tan p = 20. 


A t = 1500 GeV, 
M 2 = 200 GeV, 
M a = 1000 GeV, 
m susy = 1000 G^V , 

( 8 ) 


In addition we have the parameter m t , the soft SUSY 
breaking parameter of the top squark, which is varied in 
Fig. [6] The default value m t = 400 GeV in combination 
with S0FTSUSY [ 351 leads to the following values for 
the DR masses 


m~ h = 370 GeV, m~ h = 1045 GeV, 
m~ q = 1042 GeV, m~ g = 860 GeV, (9) 

where m q corresponds to the average of m s , mj, m~ s , m 5 
and m b and the renormalization scale has been set to the 
on-shell top quark mass. 

In Ref. 11201 the program H3m [20 211 has been 
used in order to compute the lightest MSSM Higgs 
boson mass. Combining H3m with version 2.6.5 of 
FeynHiggs m and version 3.1.1 of S0FTSUSY j35j 
leads to a Higgs boson mass of approximately 126 GeV 
almost independent of m t HI201 . 

In Fig. pla) the quantity oy QCD is shown as a func¬ 
tion of m t at FO, NFO and NNFO (from bottom to top). 
For each order three curves are shown where the dotted 
curve corresponds to the SM. The SUSY QCD correc¬ 
tions are included in the dashed and solid line where for 
the former the soft and hard renormalization scales, ji s 
and jih have been identified with M b l2 and for the latter 
p s = Mhl 2 and p b = M t has been chosen. 

One observes that the difference between SM and 
MSSM becomes small for increasing m t which is ex¬ 
pected since in this limit the spectrum becomes heavy. 
However, for smaller values of m t a sizeable effect of 
the generic supersymmetric contribution is visible. For 
example, for m t = 400 GeV a reduction of the SM cross 
section of about 5% is observed when including NNFO 
supersymmetric corrections. 

The difference between the dashed and solid line 
in Fig. [6ja) quantifies the effect of the resummation 


of In(Af^/A/^ ) where A/heavy is a heavy mass scale 
present in the calculation of C \. It is negligible for large 
m t , however, for smaller values it can lead to a visible 
effect. 

It is interesting to note that supersymmetric three- 
loop corrections to C\ computed in Refs. II1201 [121 1 
provide an important contribution to the difference of 
the solid and dotted curve in Fig. [6ja). In fact, if 
we choose fh t = 400 GeV and identify the three-loop 
matching coefficient in Eq. Q with the SM one a re¬ 
duction of the cross section of only 3% and not 5% is 
observed. 

Fet us finally present results for cr(pp —> h + X) 
which include in addition bottom quark contributions 
up to NFO and furthermore also electroweak correc¬ 
tions. In Fig. |6jb) we show the dependence on m t at 
FO, NFO and NNFO (from bottom to top). The dot¬ 
ted curves in Fig. [6jb) correspond to the solid ones of 
Fig |6ja), i.e. they only include the top-sector contribu¬ 
tion. The inclusion of the bottom quark effects at NFO 
(cf. Eq. ([7|) leads to a reduction of about 5% as shown 
by the dashed curves. The reduction is basically inde¬ 
pendent of m t and tan/3^] Thus, even for tan/3 = 20 
the bottom quark effects are small for the considered 
scenarios and, hence, at NNFO the approximation of 
vanishing bottom Yukawa coupling is justified. The re¬ 
duction due to bottom quark effects is to a large extent 
compensated by the electroweak corrections taken into 
account multiplicatively as can be seen by the solid line 
which includes all contributions of Eq. 0. 

To conclude this subsection let us remark that in the 
recent years considerable progress in the computation 
of higher order supersymmetric corrections to the Higgs 
boson production cross section has been achieved. The 
supersymmetric NNFO corrections can affect the pro¬ 
duction cross section by a few percent in case there is a 
splitting in the top squark masses by a few hundred GeV 
and the overall scale of the spectrum is not too heavy. 
Such effects are certainly relevant once the experimental 
precision for the cross section has been reduced, in par¬ 
ticular, once there are hints for new particles from direct 
searches. 

4. Renormormalization group functions in the SM 

to three loops 

Renormalization group functions are fundamental 
quantities of each quantum field theory. In general, beta 
functions provide insights in the energy dependence of 


4 The dependence on tan (3 is studied in Ref. mo). 
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Figure 6: Cross section as a function of the singlet soft SUSY breaking parameter of the right-handed top squark, fh t . (a) top quark/top squark 
contribution cr t . (b) complete contribution including also bottom quark and electroweak effects as described below Eq. J7J. Figure taken from 
Ref. fl20l . 


cross sections, hints to phase transitions and can provide 
evidence to the energy range in which a particular the¬ 
ory is valid. In the recent years the beta functions of all 
SM couplings have been extended to three loops. This 
was partly triggered by the discovery of a Higgs boson 
at the LHC i ll 26 , 1271 . A precise running of the Higgs 
boson self coupling from low energies to energies of the 
order of the grand unification scale (GUT scale) is man¬ 
datory in order to make firm statements about the sta¬ 
bility of the Higgs boson vacuum. A further motivation 
is connected to the possibility of gauge and/or Yukawa 
coupling unification at high energies. Again, the renor¬ 
malization group functions are needed in order to trans¬ 
fer the knowledge of the couplings at the electroweak 
scale to the high energy scales where unification is ex¬ 
pected. 

Within the CRC/TR 9 several important contributions 
have been achieved which are summarized in the fol¬ 


lowing subsections. In particular, in Subsections 4.1 


I4.3l we describe the calculation of the SM beta functions 
and review the important contributions to the topic. In 
Subsection |4.4| the calculation of the three-loop SUSY 
QCD beta function is mentioned. In Sections |4.5f - 
|4. 7 1 decoupling constants, which establish the relations 
between couplings and parameters in the full and ef¬ 
fective theories, are discussed. Furthermore, three ap¬ 
plications, both of the decoupling constants and the 
renormalization group functions, are presented in Sec¬ 
tions |478j [479]and[47TO| Further details can also be found 
in Ref. (22l . 


4.1. SM couplings and definition of the (3 functions 
The SM is a product of the three subgroups Uy(l), 
SU(2 ) l and SU(3)c. To each one a coupling constant is 


assigned which is usually denoted by g i, g 2 and g 3 , the 
gauge couplings. Often it is convenient to introduce 

g 2 

i= 1,2,3, (10) 

An 

which will be used below in the expressions for the 
beta functions. oq, a 2 and a 3 obey the following all¬ 
order relations to the fine structure constant ctqed, the 
weak mixing angle 6 W and the strong coupling constant, 
which are usually used in the SM 

5 (Tqed 

= o - 2~D ’ 

3 cos 2 0 W 

<*QED 

a 2 = —z — , 

sin 6w 

a 3 = a s . ( 11 ) 


Note that in Eq. ( pTj ) SU(5) normalization has been used 
which leads to the factor 5/3 in the equation for oq. 
Equation ( p~T] ) can as well be considered as a definition 
for uqed and 6 W . 

For each massive fermion there exists a Yukawa 
coupling to the Higgs boson. To lowest order it is given 
by 


a x 


ttQEDtrfc 
2 sin 2 6 W M ^ 


with x = e,p,r,u,d, c , s,t,b. 


( 12 ) 

where m x and M w are the fermion and W boson mass, 
respectively. 

Finally, there is the Higgs boson self-coupling A, 
which we define via the following term in the Lagrange 
density 


-Csm = ...-(4nl){H^H) 2 + ..., (13) 
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describing the quartic Higgs boson self-interaction. H 
is the Higgs doublet field in the SM. 

Note that the Yukawa couplings of the lighter fermi¬ 
ons are phenomenologically irrelevant. For this reason 
we only keep a t9 and a T different from zero and 
define a 4 = a u as = a t>, a 6 = <Tt and a 7 = X. Some of 
the formulae presented below can easily be extended to 
the more general case in an obvious way. 

Throughout this section we adopt the modified min¬ 
imal subtraction (MS) renormalization scheme. Note 
that in this scheme the beta functions are mass inde¬ 
pendent which allows us to perform the calculations in 
the unbroken phase of the SM where all particles are 
massless. 

We define the beta functions as 

1 u 2 Aj—=Pi({aj},e), (14) 

d/i z 71 


where e = (4 - d)/2 is the regulator of Dimensional 
Regularization with d being the space-time dimension 
used for the evaluation of the momentum integrals. The 
dependence of the couplings a t on the renormalization 
scale is suppressed in the above equation. 

The beta functions are obtained by calculating the 
renormalization constants relating bare and renormal¬ 
ized couplings which we define vi^] 

ctT = H 2 e Z ai ({aj}, e) ai . (15) 

where i = 1,..., 7. Taking into account that a ^ does 
not depend on yu, Eqs. ( [14] ) and fj"5] ) lead to 



( 16 ) 


where the first term in the first factor of Eq. © ori¬ 
ginates from the term /i 2e in Eq. (15) and vanishes in 
four space-time dimensions. The second term in the 
first factor contains the beta functions of the remaining 
six couplings of the SM. 

From the second factor of Eq. ( fib] ) it is obvious that 
three-loop corrections to Z a . are required for the com¬ 
putation of Pi to the same loop order. 

Note that for the gauge couplings a \, a 7 and the 
one-loop term of Z a . only contains or*, whereas at two 
loops all couplings are present, except X. The latter ap¬ 
pears for the first time at three-loop level. As a con¬ 
sequence, for the three-loop calculation of P\, P 2 and 


5 Note that, in the case of a 7 = A we have that Zaj contains terms 
proportional to 1 /A. The developed formalism, in particular Eq. (Th) 
is nevertheless applicable. 


P 3 , it is necessary to know Pj for j = 4,5,6 to one- 
loop order and only the £-dependent term for ^ 7 , namely 
Pi = -ea 7 / 7 r. 

This is different for the Yukawa couplings. Here, 
in the one-loop corrections to the renormalization con¬ 
stants au • • •, &6 appear which means that the two-loop 
gauge coupling beta functions are needed in order to 
compute the three-loop term to the Yukawa beta func¬ 
tion. X is present in the Yukawa coupling renormaliza¬ 
tion constant starting from two loops and thus the one- 
loop term of p 7 is needed. 

For the calculation of p 7 = p% all beta functions ex¬ 
cept /?3 are required to two-loop accuracy since the cor¬ 
responding couplings are already present in the one- 
loop renormalization constant. The strong coupling 
only enters at two loops and, hence, for ^3 only the one- 
loop result is needed. 

Before presenting some details on the calculation of 
three-loop corrections to the renormalization constants 
in the next subsection we want to end this subsection 
with a summary of important milestones for the calcu¬ 
lation of the beta functions in the SM: 

• The one-loop beta functions in gauge theories 
along with the discovery of asymptotic freedom 
have been presented in Refs. H12811129 1. 

• The corresponding two-loop corrections 

- in gauge theories without fermions 1 1301 

USE 

- in gauge theories with fermions neglecting 
Yukawa couplings III321413411 , 

- involving Yukawa corrections OSl , 
are also available. 

• The two-loop gauge coupling beta functions in 
an arbitrary quantum field theory have been con¬ 
sidered in Refs. 1 1361137 1. 

• Two-loop corrections to the renormalization group 
functions for the Yukawa and Higgs boson self¬ 
couplings in the SM are also known H137fll42l . 

• The contribution of the scalar self-interaction at 
three-loop order has been computed in 11431 1441. 

• The gauge coupling beta function in QCD to three 
loops is known from Refs. 197. 98|. 

• The three-loop corrections to the gauge coupling 
beta function involving two strong and one top 
quark Yukawa coupling have been computed in 
Ref. nn. 
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• The three-loop corrections for a general quantum 
field theory based on a single gauge group have 
been computed in 11461 . 

• The complete three-loop corrections to the gauge 
coupling beta functions in the SM has been com¬ 
puted in Refs. 1147U1491 

• The four-loop corrections in QCD are known from 
Refs. fBOlfBn . 

• The dominant three-loop corrections to the renor¬ 
malization group functions of the top quark 
Yukawa and the Higgs boson self-coupling have 
been computed in Ref. C52). In that calculation 
the gauge couplings and all the Yukawa couplings 
except the one of the top quark have been set to 
zero. 

• Complete three-loop corrections to the Yukawa 
coupling beta functions in the SM has been com¬ 
puted in Ref. Il53j (see also Ref. Il54l ). 

• Complete three-loop corrections to the Higgs bo¬ 
son self-coupling beta function in the SM has been 
computed in Refs. H55l 11561 . even with complex 
Yukawa matrices E3. 


4.2. Calculation of the renormalization constants to 
three loops 

In order to compute the renormalization constant of a 
coupling one has to consider loop corrections to a ver¬ 
tex involving this coupling. In addition the wave func¬ 
tion renormalization constants for the external particles 
have to be computed, see, e.g., Refs. [42, 145 1. For ex¬ 
ample, if we consider the Appoint vertex with external 
fields 0i,..., f n and denote its coupling constant by g , 
we obtain 


Z, = 


^01 '"(pN 




'01 ■ ■ ■ ^(pN 


(17) 


where the Z (f)i are the wave function renormalization 
constants for the 0/, Z^...^ is the corresponding ver¬ 
tex renormalization constant, and Z g the renormaliza¬ 
tion constant for the coupling g. Formulae can be de¬ 
rived where the MS Z factors are obtained from the 
ultraviolet-divergent part of amputated Green’s func¬ 
tions (accompanied by higher order e terms of lower- 
order contributions; see, e.g., Ref. EEl). 

Sample Feynman diagrams for the case of the gauge 
couplings are shown in Fig. [7] The renormalization 
constants for the couplings g\ and g 2 can, e.g., be 
computed from the gauge boson two-point function 


~~C>—0—©—©—©~~ 


-PTv -C\~ 

V.U bu*i 







X1»3>X>IK>- 


x>~ >‘-9- ]xj} TT> W 


Figure 7: Sample Feynman diagrams contributing to the Green’s 
functions which can be used for the calculation of the gauge coupling 
renormalization constants. Solid, dashed, dotted, curly and wavy lines 
correspond to fermions, Higgs bosons, ghosts, gluons and electroweak 
gauge bosons, respectively. 


(first line), ghost two-point function (third line) and the 
gauge boson-ghost vertex (fourth line). In the case of 
the Yukawa coupling the Higgs boson-fermion vertex 
can be used (together with the corresponding two-point 
functions) and for the Higgs boson self coupling the ver¬ 
tices involving four scalar particles. 

In the practical evaluation of the loop integrals one 
exploits the fact that in the MS scheme the renormaliza¬ 
tion constants do not depend on the kinematical quantit¬ 
ies like external momenta or particle masses. Thus, it is 
possible to choose convenient configurations which al¬ 
low for a simple evaluation of the integrals. In case two- 
or three-point Green’s functions have to be computed it 
is convenient to treat all involved particles as massless 
and keep only one external momentum non-zero. This 
leads to massless propagator-type integrals which up to 
three loops can be computed with the help of the FORM 
program MINCER EH). This approach can in principle 
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rapj—+ MAIAD/MINCER 


T 


FeynAits model file 


T 


FeynRulea 


Figure 8: Overview of our automated setup. Calling up the programs 
in the uppermost line determines and evaluates a given process in a 
given model. The vertical workflow leads to the implementation of a 
new model in the setup. 


lead to infrared divergences. However, introducing a 
small mass as potential infrared regulator in combina¬ 
tion with asymptotic expansion it is straightforward to 
check infrared safety, see Ref. 11481 . This method has 
been used to compute the gauge and Yukawa coupling 
renormalization constants up to three loops. 

In case the renormalization constants have to be ex¬ 
tracted from four-point functions, like the renormaliza¬ 
tion constant of the Higgs boson self coupling, the de¬ 
scribed procedure cannot be applied since nullifying all 
but one external momenta inevitably leads to infrared 
divergences. Thus, it is more profitable to set all ex¬ 
ternal momenta to zero and introduce a common mass 
M to all particles. Up to three loops the resulting loop 
integrals are well studied in the literature and a auto¬ 
mated calculation is possible with the help of the FORM 
program MAT AD EH- A minor disadvantage of this ap¬ 
proach is that additional non-standard counterterm con¬ 
tributions have to be introduced for the mass parameter 
M. A details description of this method can be found 
in in Ref. *1591 . It has been used to compute three-loop 
corrections to the anomalous dimension matrix neces¬ 
sary for analyzing the decay B —> X s y at NNLO 11601 . 

In general a large number of Feynman diagrams is 
involved in the calculation of three-loop SM Green’s 
functions ranging up to a few millions for the Higgs bo¬ 
son four point functions. Thus, an automated setup is 
mandatory. 

In Refs. 11471 148 ] a well-tested chain of programs 
has been used that work hand-in-hand: QGRAF cm 
generates all contributing Feynman diagrams. The 
output is passed via q2e 11621 11631 . which trans¬ 
forms Feynman diagrams into Feynman amplitudes, to 
exp 1162111631 that generates FORM 11641 code. The 
latter is processed by MINCER Ifl58l and/or MATAD d 
that compute the Feynman integrals and output the 


6 expansion of the result. The parallelization of 
the latter part is straightforward as the evaluation of 
each Feynman diagram corresponds to an independ¬ 
ent calculation. The input for QGRAF and q2e is 
provided by the program FeynArtsToQ2E which trans¬ 
lates FeynArts *1551 model files into model files pro- 
cessable by QGRAF and q2e. Furthermore, it is possible 
to apply the package FeynRules 1166] in order to gen¬ 
erate model files for FeynArts. The complete workflow 
is illustrated in Fig. [8] 

A similar level of automation has been obtained in 
Refs. 1149111531 where LanHEP 11671 . FeynArts DHL 
MINCER *T58l . color lfl68l and DIANA lfl69l has been 
used. 

In Ref. Q55) also QGRAF has been used for the dia¬ 
gram generation. The output is further processed with 
GEFICOM 1170M1721 and the resulting three-loop integ¬ 
rals are computed with MINCER and MATAD^] 

For some of the renormalization constants, in partic¬ 
ular the ones related to the Yukawa couplings, a careful 
treatment of the y$ matrix is required. A practical pre¬ 
scription based on the work 1 1741 for the computation 
of renormalization constants is given in Ref. 11751 and 
has been adopted in Refs. 1149, 1531 1551. 

4.3. Numerical results 

We refrain from displaying analytic results for the 
beta functions which can be found in the original public¬ 
ations. Rather we briefly discuss the numerical impact 
of the three-loop term. In Fig.[9]the running of the three 
gauge couplings is shown assuming that the SM is valid 
up to high scales. In panel (a) the energy varies over 16 
orders of magnitude. It is interesting to have a closer 
look to the intersection region of oq and oq which is 
shown in panel (b) where the one-, two- and three-loop 
results are shown as dotted, dashed and solid lines. The 
bands around the three-loop curves reflect the numerical 
uncertainties for oq(M z ) and aq(4^z) which are given 

bjQ 

afS(Mz) = 0.0169225 ±0.0000039, 

afS(Mz) = 0.033735 ±0.000020, (18) 

Defining the difference between successive loop or¬ 
ders as remaining theoretical uncertainty one observes 
that without the three-loop corrections the theory uncer¬ 
tainty is much larger than the experimental one. How¬ 
ever, after including the three-loop term the experi¬ 
mental error dominates over the theoretical one. 


6 Some of the program packages, which are not publicly available, 
can be obtained from the authors upon request; see also fl73l . 

7 See Ref. 11481 for more details. 
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Figure 9: (a) The running of the gauge couplings at three loops. 

The curve with the smallest initial value corresponds to a\ (blue), the 
middle curve to a ?2 (green), and the curve with the highest initial value 
to cr 3 (red), (b) Magnification of the intersection region of a\ and 
qt 2 where the dotted, dashed and solid lines correspond to one-, two- 
and three-loop precision, respectively. The bands around the three- 
loop curves visualize the experimental uncertainty. (Figure taken from 
Ref. fl48l .) 


As a further example we show in Fig. 10 the running 
of the quartic coupling up to the Planck scale with ini¬ 
tial conditions taken for p = M t . One observes quite 
significant effects when going from one- to two-loop 
accuracy. The three-loop corrections only lead to a 
small shift indicating a stabilization of the perturbat¬ 
ive expansion. Notably, the theory uncertainty due to 
the running is now negligible as compared to the para¬ 
metric one which is dominated by the top quark mass. 
In Fig. 10 the corresponding effect is shown as dotted 
lines. From Fig. [lO] one deduces that A becomes negat¬ 
ive for ln(yu/GeV) > 10.36 1 1761 1. Further details can be 
found in Subsection 14. 101 



Figure 10: Evolution of A = \nA with one- (brown, dashed), two- 
(blue, dash-dotted) and three-loop (red, solid) accuracy using initial 
conditions for /u = M t (see Ref. 17761 for details). For comparison the 
uncertainty induced by the top quark pole mass is shown as (black) 
dotted lines. Figure taken from Ref. EE 


4.4. Beta function in the supersymmetric QCD 

Before we report on the diagrammatic calculation 
of the three-loop gauge beta function in the supersym¬ 
metric QCD, few comments on the used regularisation 
scheme are in order. 

Whereas Dimensional Regularization in combination 
with the MS scheme is the canonical choice for higher 
order calculations within QCD, it is less appropiate for 
supersymmetric theories since it explicitely breaks su¬ 
persymmetry. This can easily be understood by count¬ 
ing degrees of freedom of fermionic and spin-1 fields. 
Whereas the former has four degrees of freedom the lat¬ 
ter has d degrees, where d is the space-time dimension 
and d 4 4. Thus a modification of the theory is ne¬ 
cessary to render the two numbers equal. A convenient 
choice has been introduced in Ref. 023 , so-called Di¬ 
mensional Reduction. The essential difference between 
Dimensional Reduction and Dimensional Regulariza¬ 
tion is that the continuation from 4 to d dimensions is 
made by compactification or dimensional reduction. In 
this scheme, the momentum (or space-time) integrals 
are d-dimensional in the usual way, whereas the num¬ 
ber of field components remains unchanged equal to 
four, and consequently supersymmetry is undisturbed 
(see also Refs. 1 178 , 179].) In practice it is conveni¬ 
ent to implement Dimensional Reduction by introdu¬ 
cing a new particle, the so-called s scalar, that account 
for the additional 4 - d components of the gauge bo¬ 
son fields. In this way, the well-established rules for 
computing momentum intregrals in Dimensional Regu¬ 
larization can also be applied to calculations performed 
in Dimensional Reduction. 
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An explicit calculation of the three-loop beta func¬ 
tion for supersymmetric QCD has been performed 
in Ref. Ea confirming results which have previ¬ 
ously been available in the literature 1 1461 1801- The 
main purpose of Ref. G75) was to perform a consist¬ 
ency check of Dimensional Reduction as regularisation 
scheme for supersymmetric theories at three-loop or¬ 
der. For this goal, several renormalization constants 
in supersymmetric QCD has been computed to three- 
loop accuracy. It has been shown that the same beta 
function is obtained from all three-particle vertices in¬ 
volving gluons, gluinos and e-scalars. The results of 
Ref. |I75) explicitly demonstrates the consistency of 
Dimensional Reduction with supersymmetry and gauge 
invariance, an important pre-requisite for (high-order) 
precision calculations in supersymmetric theories, like 
the three-loop corrections to the Higgs boson mass, cf. 
Section. [2] 

As a by-product of the calculation in 1 1751 . the pre¬ 
dicted relation between the gauge beta function and 
the gluino mass anomalous dimension was verified to 
three-loop order. In addition, the three-loop results for 
the quark mass anomalous dimension 1 1811 1 were con¬ 
firmed. 

The use of Dimensional Reduction (in contrast to 
Dimensional Regularization) is one of the technical 
challenges of the calculation. A further difficulty, 
which is discussed in detail in Ref. cm, is the treat¬ 
ment of 75 which only plays a sub-leading role in 
the gauge coupling beta functions of the SM. Another 
technical complication is connected to the Majorana 
nature of the gluino which is not treated consistently in 
QGRAF 11611 and thus a program, majoranas.pl, has 
been developed which implements the prescription of 
Ref. 1182 1 and adjusts the output of QGRAF. 

Let us for completeness present the result for the su¬ 
persymmetric beta function up to three-loop order. De¬ 
fining 


oSQCD 


fe) = 


■zr 

n> 0 


ry \H + 2 

-) / 
n / 




SQCD 


(19) 


we obtain for the first three coefficients 


qSQCD 1 rji 

Po = 4 Ca ~2 T/ ’ 

?r=\c\-T f ^-c F+ \c A \, 


nSQCD 

Pi 


|c> + r;(5c f + lc. 


+ 7 >(j c H 1 C ‘ C '~T6 C ‘ 


( 20 ) 


where Cp = (n 2 c - l)/(2n c ), Ca = n c are the quadratic 
Casimir invariants for S\J(n c ), and 27/ = rif is the num¬ 
ber of quark flavors (which is equal to the number of 
squark flavors in supersymmetric QCD). 

4.5. Decoupling of heavy particles 

The anomalous dimensions considered in the previ¬ 
ous subsections depend on the active degrees of freedom 
of the theory. In QCD, which for the following general 
discussion shall be used as a sample theory, this depend¬ 
ence is simply due to ny, the number of quarks contrib¬ 
uting to the running. In general, nf changes when fla¬ 
vour thresholds are crossed while increasing or lower¬ 
ing the energy scale. Note that it is not sufficient to 
simply raise or lower ny in the coefficients of the anom¬ 
alous dimensions and continue the running of the cor¬ 
responding parameters using the new beta function, but 
a careful construction of the low-energy effective the¬ 
ory is necessary. In the case of QCD this construction 
is straightforward since the effective theory with ny — 1 
active quark flavours has the same structure, i.e., con¬ 
tains the same operators as the full theory. The fields, 
masses and couplings in the two theories are related by 
so-called decoupling constants. A detailed discussion 
and the explicit construction of the effective theory can 
be found in Refs. [41, 42]. In these references formulae 
for the decoupling constants are provided which allow 
to compute the n-loop decoupling constants from light- 
particle Green’s functions with vanishing external mo¬ 
menta. The formalism resembles the renormalization 
procedure in the MS scheme. However, in contrast to 
the renormalization constants, the decoupling constants 
also contain finite (and also higher order e) terms. 

We want to remark that the necessity of introducing 
decoupling constants originates from the use of a mass- 
independent (here MS) renormalization scheme. As a 
consequence the strong coupling constant a s is not a 
physical quantity in the sense that it is not defined via 
a Green’s function which (at least in principle) can be 
measured. Furthermore, a s (p) is not a continuous func¬ 
tion of p but has finite steps at the energy scale where 
the heavy quark is integrated out. In Refs. 1 183 , 1841 a 
physical definition of the strong coupling has been intro¬ 
duced in the so-called momentum subtraction (MOM) 
scheme and it has been shown that running in this 
scheme, where the decoupling of heavy particles is 
automatic, is equivalent to the running and decoupling 
procedure in the MS. Let us in the following briefly 
comment on the comparison of the MOM and MS coup¬ 
lings before we continue our considerations in the MS 
scheme. 
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Figure 11: 1 /a s as a function of /u. The (blue) upper solid line con¬ 
taining a step for n = 2M^ corresponds to the MS result and the (red) 
lower solid line to the result in the MOM scheme. The (pink) dashed 
line represents 1 / a s in the MOM scheme. The (black) dotted lines 
lying on top of the MOM and MOM result correspond to the results 
obtained from MS value of a s using the proper conversion formulae 
(see Ref. 11841 ). Figure taken from Ref. 11841 . 


In Fig. 11 we show that the two schemes are equival¬ 
ent by plotting the inverse stron g coup fing as a function 
of yu both for the MS, MOM and MOIV^] schemes, where 
in all cases the three-loop approximation is used for the 
running and the conversion between the schemes. We 


choose aP(M z ) as the input quantity and convert at this 
scale to the other two schemes. The evolution of the 
MS coupling to lower ju values is shown by the upper 
solid lines with a step at the value fi = 2 Mb since at 
that scale the bottom quark is decoupled. Numerically 
very close is the dashed curve representing the MOM 
scheme result. The lower solid line represents the res¬ 
ult in the MOM scheme. Both for the MOM and MOM 
results, the conversion is performed for /i = M z , and the 
running to other values of fi is achieved using the corres¬ 
ponding ft function. The dotted lines on top of the MOM 
and MOM curves represent the results where the trans¬ 
formation from the MS values is performed just at the 
considered value of fi. This shows that the MS scheme 
including the described decoupling procedure is equi¬ 
valent to a physical scheme with a physical definition of 
a s . 


Let us now return to QCD with a s defined in the MS 
scheme. The phenomenologically most important de¬ 
coupling constants are the ones of a s , g as , and the light 
quark masses, £ mq . Two-loop corrections to g as have 
been computed for the first time in Refs. 11851 186], 
however, relatively complicated integrals had to be 
solved. For example, in Ref. gsd a three-loop calcu- 


8 The MOM scheme is a version of the MOM scheme which leads 
to numerical values for a s close to the ones in the MS scheme. 


lation for the Z boson decay rate has been performed. 
In Ref. E3 the two-loop result has been checked and 
the three-loop results for g as and £ mq has been added. 
The latter required the computation of three-loop one- 
scale vacuum integrals. The four-loop corrections to the 
decoupling constants have been obtained by two inde¬ 
pendent calculations performed in Refs. [44. 45]. 

The loop-order used for the running and the one used 
for the decoupling are related. In fact, n-loop running 
requires (n - l)-loop decoupling relations. This can 
be seen by considering a physical quantity R for which 
the perturbative expansion is know up to order a n s . Of 
course, in such a situation one would apply n-loop cor¬ 
rections to the beta function. However, using n -loop 
decoupling relations would affect the t erm °f R 

which is beyond the considered loop-order. 


Let us in the following demonstrate the effect of 
higher order corrections to the running and decoup- 

/'3'v 

ling by considering the relation between ar s (M r ) and 
af\M z ) (see also Refs. ATI l44l l45l [[87) where the 
charm and bottom flavour threshold is crossed at the 
scales ji c and /Jb, respectively. Note that these scales are 
not determined by theory. On general grounds one ex¬ 
pects that the result for aP(M z ) gets more and more in¬ 
sensitive on the precise choice of the decoupling scales 
when including higher order corrections. The proced¬ 
ure to compute aP(M z ) from aP(M r ) is as follows. 
In a first step we calculate aP(pi c ) by integrating the 
beta function to n -loop order with the initial condition 
aP(M T ) = 0.332. Afterwards aP(/J c ) is obtained from 
aP(ji c ) = aP(ji c )£ as using the (n - l)-loop approx¬ 
imation for £ ag . Next n-loop running is used to obtain 
a P(f*b) and the decoupling procedure is applied in ana¬ 
logy to the charm threshold to arrive at aP(jj.b)- Finally, 
we compute aP(M z ) using again the n-loop beta func¬ 
tion. For the on-shell charm and bottom quark masses 
we use M c - 1.65 GeV and M \ - 4.75 GeV, respect¬ 
ively. These values are obtained by using three-loop re¬ 
lations between the MS and on-shell quark masses [E5T- 

na. 


In Fig.[12lthe result for aP(M z ) for fixed jn c = 3 GeV 


as a function of fib is displayed for the one- to five-loop 
analysis. For illustration, fib is varied rather extremely, 
by about two orders of magnitude. While the leading- 
order result (upper right dotted line) exhibits a strong 
logarithmic behaviour, the analysis is gradually getting 
more stable as we go to higher orders. The five-loop 
curve is almost flat for //^ > 1 GeV (note the scale on the 
y axis) and demonstrates an even more stable behaviour 
than the four-loop analysis of Ref. BTl . It should be 
noted that around ^ « 1 GeV both the three-, four- and 
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Figure 12: ^ dependence of afp(Mz) calculated from af*\M T ). The 
procedure is described in the text. The dotted, steep dashed, lower 
dashed, and dash-dotted line corresponds to one- to four-loop running. 
The solid curve and the dashed line slightly below it includes the effect 
of the four-loop matching term. 


five-loop curves show a strong variation which can be 
interpreted as a sign for the breakdown of perturbation 
theory. Besides the pt dependence of af\M z ), also its 
absolute normalization is significantly affected by the 
higher orders. At the central matching scale pb - Mb, 
we encounter a rapid convergence behaviour. 

Since the five-loop coefficient of the QCD beta func¬ 
tion is not yet known we choose the (^/-independent) 
values = 0 (solid line) and /^ n/) = 150 (dashed line 
parallel to the solid on eQ Larger values would bring 
the five-loop curves even closer to the four-loop one. 

From Fig. [l2]it is possible to estimate a theory un¬ 
certainty on wp(M z ) as obtained from af \M T ) due 
to missing higher order corrections. If we restrict 
ourselves to a range of Hb between 2 GeV and 10 GeV 
and take the difference between the three- and four- 
loop curve as an estimate for the uncertainty we obtain 
Sa ( p(M z ) ~ 0.0002. The difference between the four- 
and (dashed) five-loop curve would lead to 6af\M z ) ~ 
0.0003. The variation of af\M z ) due to the variation 
of [ib leads to an additional uncertainty of 6a^\M z ) « 
0.0002. A similar uncertainty is obtained from the vari¬ 
ation of [i c between 2 GeV and 5 GeV. (This can eas¬ 
ily be checked with the program RunDec 1118811 .) Thus 
a total uncertainty of ±0.0004 (obtained by adding the 
three uncertainties in quadrature) should be assigned to 
af\M z ). The uncertainties induced by the errors in the 
quark masses are much smaller. 


9 The normalization corresponding to {po,(3\,p2,P?>} 
{1.92,2.42,2.83,18.85} for n f = 5. 


Note that it is straightforward to reproduce Fig. [12 
using the program RunDec [ 1881 which is written in 
Mathematica or CRunDec 11891 written in C++. 

The formalism which has been described above al¬ 
lows the decoupling of one heavy quark at a time. This 
procedure is certainly justified for the top quark being 
more than factor 30 heavier than the bottom quark. On 
the other hand, the ratio between the bottom and charm 
quark mass is only approximately a factor three. For this 
reason the simultaneous decoupling of the charm and 
bottom has been studied in Ref. 11901 and a decoupling 
constant relating the strong coupling defined with three 
active flavours to the one in the five-flavour theory has 
been derived. These results can be used in order to study 
the effect of power-suppressed terms in M c /Mb which 
are neglected in the conventional approach lf4ll . Vari¬ 
ous analyses are performed which indicate that the mass 
corrections present in the one-step approach are small 
as compared to lo g(p 2 /M^ b ) which are resummed using 
the conventional two-step procedure, and, thus, even for 
the charm and bottom quark case two-step decoupling 
is preferable. 


4.6. Decoupling at the SUSY scale 

As in QCD also in its supersymmetric extension, it is 
convenient to use a mass-independent renormalization 
scheme. Thus, by construction, the beta function gov¬ 
erning the running of a s is independent of the particle 
masses. It only depends on the particle content of the 
underlying theory, i.e., the number of active quarks, 
squarks and gluinos. 

In the MSSM, the one-loop decoupling relation for 
a s has been computed in Refs. 1 109 , 1911 and the two- 
loop relation is known from Refs. 1192111941 (see also 
Ref. 11951 ). Due to many different mass scales ex¬ 
act results cannot be obtained at three loops. Thus, in 
Ref. 1121 1 three different hierarchies among the super- 
symmetric particle masses have been assumed to com¬ 
pute the decoupling relation. 

Let us in the following demonstrate the numerical ef¬ 
fect of higher order corrections to the decoupling con¬ 
stants by repeating the considerations of Subsection |4.5| 
However, instead of considering the strong coupling 
with three and five active flavours we study the relation 
between af QCD) (MouT) (defined in the DR scheme) 
with M gut = 2 • 10 16 GeV and af(M z ) (in the MS 
scheme) as a function of the decoupling scale p^ tc . At 
this scale all supersymmetric particles are integrated out 
and the transition from the MSSM to the SM is made. 
We proceed as follows: In a first step we run in the SM 
from p - M z to p - p^ QC where the decoupling of 
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Figure 13: af' ^ CD) (Mgut) as a function of ^dec- Thick and thin lines 
correspond to the one- and two-step scenario, respectively (see text). 
Thin lines are only shown for three- and four-loop running. (Figure 
taken from Ref. 11211 .1 


the top quark and the supersymmetric particles is per¬ 
formed simultaneously and a^p(ji dec) is transformed to 
a ( v SQCD) 0u dec ). The use of the SUSY QCD p function fi- 


13 


nally leads to af { QCD) (M G ut)- The thick lines in Fig. 
correspond to this procedure where dotted, dash-dotted, 
dashed and solid lines correspond to one-, two-, three- 
and four-loop running and the use of decoupling con¬ 
stants to one order less. For a complete list of input 
parameters we refer to Ref. 1112111 . 

As expected, the inclusion of higher order corrections 
reduces dramatically the dependence on /i de c leading to 
an almost flat behaviour at four loops, even though the 
decoupling scale is varied over more than two orders 
of magnitude. The band around the four-loop curve 
reflects the uncertainty of the strong coupling constant 


which has been chosen as af\M z ) = 0.1184 ± 0.0007. 
It is interesting to note that // dec = M z is in general a 
less favourable choice unless three- or four-loop run¬ 
ning is used. On the contrary, for // dec « 1 TeV, which 
approximatley corresponds to the masses of the super- 
symmetric particles, higher order corrections are quite 
small. 

Alternatively, in order to obtain the thin lines we in¬ 
tegrate out the top quark in a separate step at the scale 
p - M t (M t is the on-shell top quark mass) and trans¬ 
form afterwards af\M t ) to af^ CD \M G ut) in analogy 
to the previous description. The three- and four-loop 
curves show in this variant an even flatter behaviour. 


4. 7. Decoupling at the GUT scale 

Another type of threshold corrections that have to be 
considered in the context of GUTs are those generated 
by the super-heavy particles present in such models. In 


particular, the gauge coupling unification (one of the 
most important predictions of a GUT) is very sensit¬ 
ive to these corrections. This property enables us to 
constrain the models, once the low energy values of 
the gauge couplings and their evolution are known pre¬ 
cisely. In the following, we restrict our discussion to the 
effects of the threshold corrections at the super-heavy 
scale on the energy evolution of the SM gauge coup¬ 
lings. 

The new features of the GUT threshold corrections 
as compared with those discussed in the previous sub¬ 
sections are related to the spontaneous breaking of the 
gauge symmetry. Explicitly, they take into account the 
effects of the spontaneous breaking of the GUT gauge 
group to the SM gauge group. In consequence, the 
three SM gauge couplings are affected differently by 
these corrections, so that after passing the super-heavy 
threshold they become equal and evolve as a unique 
coupling towards the Planck scale. Furthermore, the 
calculation of the associated decoupling constants is 
much more involved and very model dependent. An¬ 
other subtle point concerns the choice of gauge within 
the full theory so that the gauge invariance of the effect¬ 
ive theory, obtained by integrating out the heavy degrees 
of freedom, is also maintained 11961 . 

Currently, the super-heavy threshold corrections to 
the gauge couplings are known at the one-loop level for 
a general model 1191,11961. However, at the two-loop 
level there are only the attempts of Refs. 11971 1981 to¬ 
wards a general derivation. Precisely, this computation 
covers only theories without super-heavy fermions and 
with a rather simplified structure for the GUT breaking 
scalars. Also, trilinear interactions in the scalar poten¬ 
tial are not considered. Thus, their applicability to spe¬ 
cific models is rather restricted. 

For exemplification, we show in Fig. pA| sample Feyn¬ 
man diagrams contributing at the two-loop level to 
the decoupling constant of the gauge couplings. As 
usual, to compute the decoupling constants one con¬ 
siders Green’s functions involving only light external 
particles and at least one heavy field within the loops. 
In Fig. [14] the contributions to the self-energies of the 
light gauge bosons and their associated ghosts, and to 
the gauge boson-ghost three-point function are shown. 

A numerical analysis of the effect of the two-loop 
GUT threshold corrections on gauge coupling unifica¬ 
tion was performed in Ref. 11971 . It turned out that, 
for models containing large representations, the effects 
of the two-loop GUT threshold corrections exceed by 
almost an order of magnitude those induced by the cur¬ 
rent experimental uncertainties of the gauge coupling 
constants. In consequence, the GUT threshold correc- 




















L.N. Mihaila, M. Steinhauser 


20 


vwyjwvw\|ww wv| 


VW\^|^VWV VW^^^^VWV 



Figure 14: Feynman diagrams contributing to the two-loop threshold 
corrections at the GUT scale. Bold (red) lines represent fields with 
mass of O(Mquj) and thin lines massless fields. Furthermore, curly 
lines denote gauge bosons, dotted lines ghosts, dashed lines scalar 
fields and solid lines fermions. 


tions necessarily have to be taken into account in any 
phenomenological study concerning the gauge coup¬ 
ling unification. An impression of their numerical ef¬ 
fects in the minimal SUSY SU(5) model can be seen in 
Fig.[l6|b) (however, only at the one-loop level). 


4.8. Application I: low-energy theorems 

There is a close connection between the decoupling 
constants and the effective coupling of a Higgs boson 
to gluons and light quarks which at first sight is quite 
surprising. It is established by the so-called low-energy 
theorem (LET) which relates the decoupling constants 
to the Wilson coefficients in the effective Lagrangian 


X e ff = |ci+ C2mb^b^bj , (21) 

where the first part has already been shown in Eq. 0- 
The term proportional to C 2 has been specified to the 
bottom quarks; similar contributions also exist for the 
other light quarks. In Ref. ED the following LETs have 
been derived 


Ci 


C 2 


m,d 

dm, 


u 


l + ^-r 
dm A h - 


( 22 ) 


where we refer to ED for a proper definition of the de¬ 
coupling constants (with g as = (£ g ) 2 ). 

At lowest order the LETs given in Eq. ( |22| ) can easily 
be motivated. Consider, e.g., C\ which can be com¬ 
puted from the one-loop top quark triangle with two 
gluons and a zero-momentum Higgs boson as external 
particle. On the other hand, only the one-loop heavy- 
top diagram of the gluon propagator contributes to g as . 
Zero-momentum top quark-Higgs boson coupling are 
obviously generated by taking derivatives with respect 
to the top quark mass which leads to the first line in 


Eq. ( [22] ) expanded to one-loop order. At higher order 
this simple picture does not work any more due to the 
fact that also other Green’s functions enter e.g., the 
ghost two-point function. 

It is interesting to note that Eq. ( [22] ) only contains log¬ 
arithmic derivatives. Thus, it is sufficient to know the 
ln(m 2 ) terms of £ as and £ mq in order to obtain C\ and C 2 
at the respective loop order. Since m t is the only dimen¬ 
sionful physical scale it actually occurs in the combin¬ 
ation In (p 2 1 m 2 t ) and thus the logarithms at order a ( f +1 j 
can be reconstructed from the order a n s with the help of 
renormalization group functions. In Refs. [44, 45 ] this 
has been exploited to obtain the five-loop result for C\ 
from the four-loop result of In beyond-SM the¬ 

ories this trick cannot be applied due to the presence of 
more than one particle mass. 

For earlier work on LETs we want to refer to 
Refs. 1 1199 , [200 1 and for the application of LETs 
in supersymmetric theories we want to mention the 
works lll0lll2Hl201l . In Ref. 11901 it has been applied 
to QCD involving more than one heavy quark. 

4.9. Application II: gauge coupling unification in mod¬ 
els based on the SU(5) group 

The quantum numbers of the SM fermions together 
with the apparent convergence of the strong and elec- 
troweak couplings at energies below the Planck scale 
point towards a unified description of the SM interac¬ 
tions. Furthermore, one of the fundamental predictions 
of a GUT is the existence of baryon and lepton number 
violating interactions which can manifest themselves at 
low energy via matter instability (for a review see for ex¬ 
ample Ref. l202i ). Though the decay of the proton has 
not been observed so far, the lower bound on the proton 
lifetime together with the low-energy values of the SM 
gauge couplings and the SM fermion masses and mix¬ 
ing provide us severe constraints on the class of viable 
GUT models. 

Gauge coupling unification is highly sensitive to the 
mass spectrum. This property allows us to probe the 
unification assumption through precision measurements 
of low-energy parameters like the gauge couplings at 
the electroweak scale or the mass spectrum and preci¬ 
sion calculations. Through such analyses, we can make 
predictions for some of the mass parameters of the mod¬ 
els that have to be compared with the constraints de¬ 
rived from the non-observation of the proton decay. In 


10 Note that the fermionic contribution of the five-loop beta func¬ 
tion, which is not yet known, enters the five-loop term of Ci as a free 
parameter. 
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addition, GUTs might predict interesting signatures at 
the LHC (see, e.g., Ref. 1203 1 for a recent search for 
heavy fermionic triplets and Ref. 12041 for searches for 
scalar leptoquarks). These additional constraints from 
the LHC can, in some cases, be sufficient to even rule 
out models 1205142091 . 

In the following, we restrict our discussion to GUTs 
based on the SU(5) gauge group. Although phenomen¬ 
ologically there are other theories which are better mo¬ 
tivated, SU(5) GUTs are most predictive. For example, 
one of the few absolute certainties about grand unific¬ 
ation today is that the original SU(5) model of Georgi 
and Glashow (GG) 12101 is ruled out. In particular, the 
failure of the minimal model can be attributed partially 
to the lack of gauge coupling unification 1211112131 . 

Let us briefly recall the reason why gauge coupling 
unification fails within the minimal GG model FH While 
oq and a 3 meet around 10 16 GeV, oq and a 2 intersect 
already at about 10 13 GeV, at odds with the bounds en¬ 
forced by the nonobservation of the proton decay. More 
precisely, model independent upper bounds on the pro¬ 
ton lifetime 1 214 1 together with the latest experimental 
data from the Super-Kamiokande observatory 12151 im¬ 
ply a conservative lower bound on the unification scale 
M g of about 10 15 - 5 GeV. Hence, the key ingredients 
for a viable unification pattern are additional particles 
charged under the SU(2) L group that delay the meeting 
of a 1 and a 2 . This is, essentially, the philosophy behind 
two recent proposals where an extra scalar representa¬ 
tion 15 H ■216U217I . or alternatively, a fermionic repres¬ 
entation 24 F [2051 218 1 are added to the field content of 
the model. In both cases, the extra degrees of freedom 
have the correct quantum numbers to restore unification 
by properly modifying the running of the gauge coup¬ 
lings. 

In particular, the study performed in Ref. 12091 
proved that a three-loop analysis of the gauge coupling 
unification is required within the SU(5)+24/7 model. As 
can be read from the Fig. [15] the two-loop corrections to 
the mass of the electroweak triplets m 3 for a fixed uni¬ 
fication scale M g is of the same order of magnitude as 
the one-loop contributions and amount to several TeV. 
The three-loop corrections are rather small (hundreds of 
GeV) and lie within the uncertainty band of the two- 
loop results. Thus, the three-loop order analysis is re¬ 
quired in order to reduce the theoretical uncertainties at 
the level of parametric uncertainties induced by the low- 
energy measurements of the electroweak gauge coup¬ 
lings. Moreover, the achieved theoretical accuracy to¬ 
gether with experimental data in the reach of the LHC 

11 For the definition of a\, and «3 see Eq. 0 



log 10 (M c /GeV) 


Figure 15: The maximal value of the effective triplet mass m 3 
in the SU(5)+24/r model as a function of the unification scale Mq. 
The black, blue and red bands (from bottom-left to top-right) cor¬ 
respond respectively to the one-, two- and three-loop running ana¬ 
lysis. The error bands are due to the lcr uncertainties on the low- 
energy electroweak couplings a\ ( Mz ) and aiiMz) (Figure taken from 

Ref. (255)). 


and the Super-Kamiokande observatory will provide us 
with sufficient information to even disprove the model. 

Another example of a predictive GUT model is the 
minimal SUSY SU(5) GT9] [220) . It has the important 
feature that one can derive unambiguous correlations 
among its parameters and even rule out the model once 
sufficiently precise experimental data become available. 
Immediately after the formulation of the minimal SUSY 
SU(5) it has been noticed that within SUSY GUTs 
new dimension-five operators cause a rapid proton de¬ 
cay 11221112221 . This aspect was intensively studied over 
the last thirty years with the extreme conclusions of 
Refs. 12231^241 that the minimal SUSY SU(5) model 
is ruled out by the combined constraints from proton 
decay and gauge coupling unification. However, care¬ 
ful analyses have shown that taking into account fer¬ 
mion mixing 12251 or higher dimensional operators in¬ 
duced at the Planck scale (226, 227], one can substan¬ 
tially weaken the constraints and the minimal model 
remains a valid theory. Also the high-precision ana¬ 
lysis of Ref. 1228] confirmed that the minimal SUSY 
SU(5) model cannot be excluded by the current ex¬ 
perimental data. In particular, one observes an in¬ 
crease of the super-heavy Higgs triplet mass by about 
an order of magnitude when three-loop effects are con¬ 
sidered. These results attenuate substantially the tension 
between the theoretical predictions and the constraints 
derived from the experimental data. 

Beyond the minimal version of the SUSY SU(5) 
model, the interplay between the theoretical predic- 
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tions and the experimental data is not completely de¬ 
termined. For example, the most popular extension of 
the minimal SUSY SU(5) model, the Missing Doublet 
Model 1229112301 . designed to avoid unnatural doublet- 
triplet splitting, cannot be excluded using only the cur¬ 
rently available theoretical and experimental data. The 
model contains additional free parameters as compared 
to the minimal model and, consequently, is affected by 
large theoretical uncertainties, so that no firm conclu¬ 
sion can be drawn. 



(a) 



Figure 16: Running of the gauge couplings from the electroweak to 
the Planck scale from Ref. (228). The discontinuity for // - ususy = 
10 3 GeV and // = //gut = 10 15 5 GeV are clearly visible. In panel 
(b) an enlargement of (a) for the region around // = /tout is shown 
where for the decoupling the two values //gut = 10 16 GeV « 3.2 • 
10 15 GeV and //gut = 10 16 GeV have been chosen. Figure taken 
from Ref. [ 2281 . 


For illustration, we briefly discuss the analysis of the 
gauge coupling unification in the SUSY SU(5) model at 
three-loop level. Crucial input parameters for this ana¬ 
lysis are the precise values of the gauge couplings at the 


electroweak scale. For their determination from the ex¬ 
perimentally measured observables, one has to take into 
consideration threshold corrections at the Z-boson mass 
and at the top quark mass (for a detailed description 
see Ref. [1228 j). Furthermore, one applies the “running 
and decoupling” approach described in the previous sec¬ 
tions. The dependence on the heavy particle masses be¬ 
comes explicit through the decoupling constants. The 
constraint of gauge coupling unification translates into 
restrictions (usually expressed as correlations) on the 
mass spectrum. 

In Fig. [16] the evolution of the SM gauge couplings 
from the electroweak up to the Planck scale is shown in 
the minimal SUSY SU(5) model. The specific threshold 
corrections for this model are those due to supersym¬ 
metric particles at the TeV scale and those due to super¬ 
heavy particles at around 10 16 GeV. The associated de¬ 
coupling scales have been chosen at psusy = 1000 GeV 
and //gut = 10 16 GeV. One can clearly see the dis¬ 
continuities at the matching scales and the change of 
the slopes when passing them. The (blue) band around 
the strong coupling corresponds to the present exper¬ 
imental uncertainty. In panel (b) the region around 
ji - 10 16 GeV is enlarged which allows for a closer 
look at the unification region. Here, the decoupling of 
the super-heavy particles is performed at two different 
values of //gut- One observes quite different threshold 
corrections leading to a nice agreement of the (common) 
gauge coupling above 10 16 GeV. On the one hand, the 
plot proves that the GUT threshold corrections are es¬ 
sential for a successful unification. On the other hand, 
the dependence on //gut is interpreted as a measure of 
theoretical uncertainties because this parameter is not 
fixed by the theory. Thus, the three-loop analysis is suf¬ 
ficient to cope with the current experimental precision. 

To summarise, threshold corrections at the GUT scale 
are indispensable ingredients for precision analyses of 
gauge coupling unification. In particular, they have 
the virtue to provide us with the necessary information 
about the specific model dependence. In this way, we 
are able to perform consistency tests of various GUTs 
and, in some cases, even to refute models. 


4.10. Application III: stability of the SM vacuum 
With the discovery of a Higgs boson at the LHC, the 
question of the SM vacuum stability received a renewed 
attention. This is because the mass of the Higgs bo¬ 
son in the SM is an important indicator for the pres¬ 
ence of new physics at high energy scales. It is well 
known that, if the Higg boson mass satisfies the condi¬ 
tion M meta < Mh < M Landau then the SM is a consist¬ 
ent theory from the electroweak scale up to the Planck 
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scale. The upper limit originates from the requirement 
that the Higgs self-coupling remains in the perturbat¬ 
ive regime up to the Planck scale. Explicitly, the oc¬ 
currence of the Landau pole is avoided. The estimated 
value for M Landau is around 175 GeV II231L which is 
excluded by the direct searches at the LHC and the Tev- 
atron. The lower limit is derived from instability con¬ 
straints. The existence of a Higgs boson with a mass 
smaller than M meta would imply that new physics be¬ 
low the Planck scale is required to stabilize the SM va¬ 
cuum. The numerical value for M meta was estimated to 
111 GeV 1232]. In consequence, for a Higgs boson with 
mass Mh ~ 125 GeV one distinguishes two situations: 
either Mh > M stablllty and the electroweak vacuum is 
absolutely stable, or M meta < Mh < M stablllty , which 
corresponds to a metastable vaccum with a life-time ex¬ 
ceeding that of the Universe. 

The precise determination of M stablllty was the sub¬ 
ject of numerous recent analyses l26ll233H235l . The 
absolute stability bound on the Higgs mass M st;iblllty is 
defined as the value for which the effective potential at 
the electroweak minimum, 0 e w, and at a second min¬ 
imum at large field values, 0 > 0 ew , are the same. The 
equations to be solved read 


V(^ ew ,M stability ) 

dV 

^ 0ew 


V(0,M stability ), 

dV 



(23) 


Although the effective potential as well as the positions 
of its extrema are both gauge dependent, the solution for 
^stability j s g au g e independent as shown in Ref. 12361 . 

The most advanced recent works use for the cal¬ 
culation of M stablllty the two-loop SM effective poten¬ 
tial l ll41 , 237 1 and the three-loop RGEs for the SM 
couplings 1114711149111521115311155111561 , necessary to 
resum the large logarithms that occur in the effective 
potential V for large field values. In addition, two-loop 
corrections to the relations between the MS parameters 
and the physical observables have been considered in 
Refs. 126, 1233 ,, 2351 . The phenomenological implica¬ 
tions of the new determination of M stablllty are summar¬ 
ised in the phase diagram shown in Lig. 17 The regions 
of stability, metastability and instability of the SM va¬ 
cuum are shown for Higgs boson (Mh) and top quark 
(M t ) masses in the range corresponding to the meas¬ 
ured values. Also shown are contour lines indicating 
the instability scale Asm, defined as the scale where new 
physics is required to stabilize the SM effective poten¬ 
tial. It appears that the measured values for Mh and M t 
are rather special, in the sense that, the SM vacuum lies 
at the border between stability and metastability. One 
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Figure 17: SM phase diagram in terms of the Higgs boson Mh and top 
quark M t masses. The dotted contour-lines show the instability scale 
Asm in GeV for a s (Mz) = 0.1184. Figure taken from Ref. 12351 . 




Figure 18: Instability scale as a function of £ = £w(M t ) = for 

the Fermi gauge. The dashed curve corresponds to the case where the 
gauge fixing parameters do not run. 


also concludes that a metastable electroweak vacuum 
can comply with the data and new physics below the 
Planck scale is not necessarily implied. 

Note that the interpretation of the SM instability scale 
as the scale A, where 


V((p) > V(0 ew ) for all 0 < A (24) 

is problematic. The effective potential is gauge depend¬ 
ent, and it renders A as solution of the inequality (24) 
gauge dependent as well. The numerical effects due to 
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gauge dependence are rather large as has been explicitly 
shown in Refs. I ;236 . ' 2381 [2391. The instability scale 
can vary by few orders of magnitude in 7^-like gauges, 
depending on the choice of the gauge fixing paramet¬ 
ers. This behaviour is illustrated in Fig. [18] taken from 
Ref. 12361 . Here, the dependence of the SM instability 
scale [defined as in inequality ( [24] )] on a common gauge 
fixing parameter in the electroweak sector is displayed 
for the special choice of a Fermi gauge (see Ref. [ 236 ] 
for precise definition). As can be read from the figure, 
the resummation obtained through the use of running 
gauge fixing parameters is essential to reduce the vari¬ 
ation of the instability scale. However, even the RGE 
improved prediction for Asm varies by an order of mag¬ 
nitude for gauge fixing parameters in the perturbative 
regime. Let us at this point stress that significant vari¬ 
ations of the instability scale are also observed when 
changing the form of the gauge fixing term in the Lag- 
range density. 

To conclude, a metastable electroweak vacuum can 
comply with the present experimental data, and new 
physics below the Planck scale is not necessarily re¬ 
quired. Furthermore, the physical threshold for new 
physics can be determined only after specifying an ul¬ 
traviolet completion of the SM. Its determination from 
the stability requirement of the SM potential is theoret¬ 
ically inconsistent due to the gauge dependence as dis¬ 
cussed above. 

5. Summary 

With the discovery of a Higgs boson particle at the 
LHC physics has entered a new era which is at the mo¬ 
ment dominated by detailed and careful studies of its 
properties. The central question about the particle dis¬ 
covered at about 125 GeV is whether this is “the Higgs 
boson” of the SM or only one degree of freedom of a 
bigger theory. The fundamental difficulty to identify 
new physics signatures at the LHC, so far, initiated nu¬ 
merous precision analyses focused on the Higgs sector. 

In this review we report on recent precision calcu¬ 
lations and computational developments. In particular, 
the following topics are highlighted: the determination 
of the Higgs boson mass within the MSSM at three-loop 
accuracy; the calculation of the Higgs production in the 
SM and the MSSM through N 3 LO and NNLO, respect¬ 
ively; the computation of renormalization group func¬ 
tions within the SM and its extensions at the three-loop 
order. Special emphasis was put on the realisation of the 
gauge coupling unification within BSM theories and on 
the study of the SM vacuum stability. 


Hopefully, accurate understanding of the Higgs phe¬ 
nomenology obtained through high precision analyses 
together with new information from experiments at the 
LHC will provide us a tool for exploring new physics 
and to find explanations for some of the long-standing 
questions in particle physics. 
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